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ABSTRACT 


An  expansion  of  the  density  matrix  is  given  into  irreducible 
SU(2)  tensors,  i.e.,  into  quantities  of  good  angular  momentum. 

These  irreducible  tensors  can  be  handled  by  all  the  powerful  tools 
developed  in  the  context  of  the  handling  of  angular  momentum.  As 
examples,  the  density  matrix  of  a cryogeni cal ly  alligned  nucleus 
is  derived  and  the  construction  of  the  angular  distributions  of 
nuclear  reactions  in  terms  of  density  matrices  is  demonstrated. 

Key  words:  angular  distributions;  angular  momentum;  impure  quantum 

states;  polarization;  quantum  theory;  reaction  theory 


I.  INTRODUCTION 


As  is  well  known  one  may  express  all  of  quantum  mechanics  in  terms  of 
density  matrices.  One  must  use  density  matrices  when  dealing  with  impure 
states  [ 1] . 1 Most  prominent  among  such  cases  is  the  description  of  systems  in 
thermal  equilibrium  with  some  heat  bath  at  non-zero  temperature.  Examples  for 
this  are  cryogenically  aligned  or  polarized  targets  for  nuclear  physics  exper- 
iments in  which  angular  distributions  or  polarization,  etc.,  are  studied.  As 
has  been  variously  emphasized,  such  situations  are  best  described  in  terms  of 
irreducible  tensors,  i.e.,  of  quantities  which  undergo  certain  well-defined 
transformations  under  rotations  [2].  The  aim  of  this  note  is  to  propose 
conventions  for  defining  irreducible  density  matrices  which  are  designed  to 
allow  the  incorporation  of  these  quantities  into  the  angular  momentum  calculus 
in  as  smooth  a manner  as  possible,  i.e.,  in  such  a way  that  no  superfluous 
phases  appear.  In  other  words,  this  note  aims  at  extending  an  earlier  paper 
[3]  to  include  density  matrices  into  a consistent  set  of  definitions.  This 
way  they  then  can  be  treated  by  all  the  powerful  tools  of  Racah  calculus 
including  the  graphical  method  introduced  in  Ref.  3 (see  also  Ref.  4).  We 
shall  employ  these  techniques  in  the  present  note,  and  refer  the  reader  for  a 
complete  description  of  the  notation  and  the  techniques  to  Chapter  4 of 
Ref.  5. 

In  section  II  the  density  matrix  is  introduced  in  the  context  of  the 
evaluation  of  the  expectation  value  of  an  operator  for  the  case  of  a pure 
state.  The  definition  of  the  irreducible  density  matrix  this  way  is  achieved 
in  a natural  way.  As  examples  for  the  use  of  density  matrices  we  describe  the 


lumbers  in  brackets  indicate  the  literature  references  at  the  end  of  this 
document . 
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construction  of  the  expressions  for  the  angular  distributions  and  polariza- 
tions of  reactions  of  arbitrary  complexity.  It  is  here  that  the  advantage  of 
using  the  consistent  definitions  for  the  density  matrices  becomes  apparent  in 
that  one  automatically  achieves  the  complete  expressions  including  the  inter- 
ference terms  in  a manner  which  essentially  eliminates  the  possibility  of 
errors,  including  errors  in  phasing.  The -expressions  are  written  in  terms  of 
the  density  matrices  of  the  initial  and  final  states  and  of  the  analyzing 
instruments,  and  thus  are  fully  general.  In  this  section  also  the  special 
precautions  are  mentioned  one  must  observe  if  one  works  in  the  occupation 
number  representation,  i.e.,  if  the  amplitudes  are  creation-annihilation 
operators . 

The  general  case  of  an  impure  state  is  considered  in  section  III  where  it 
is  shown  that  all  the  results  of  section  II,  i.e.,  of  the  density-matrix 
formulation  for  pure  states,  are  immediately  applicable  in  the  general  case. 

As  an  example  of  an  impure  state  the  density  matrix  of  a nucleus  in  thermal 
equilibrium  with  a crystal  lattice  is  derived  in  full  detail.  This  is  done 
for  the  case  where  the  lattice  field  is  of  arbitrary  multipolarity. 

The  relation  of  the  irreducible  density  matrix  to  the  density  matrix 
defined  in  the  usual  way  is  demonstrated  in  section  IV.  Finally,  in  the 
Appendix,  some  definitions  are  collected  for  the  convenience  of  the  reader. 
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II.  PURE  STATES 


Consider  a state  of  given  angular  momentum  and  polarization.  It  is 
written  as  the  invariant 


l 

M 


,(» 


I =:  (21  + 1)* 


(1) 


We  use  the  usual  convention  [2-5] 


U)  =.  (-)I+M  a[I] 


‘M 


-M 


(2) 


Care  has  to  be  taken  when  constructing  the  "complex  conjugate"  state  if  it 
describes  a system  containing  Fermions.  We  shall  therefore  separate  the 
conjugation  explicitly  into  the  two  components:  change  of  the  gender 
(contrastandard  — costandard)  and  "transpose."  We  thus  write 


f -l  a.(I)f 

M 


M 


II]t 


(3) 


together  with  the  definitions  (convention): 


r*  ■=  4n  = (-)I+M^] 
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(4a) 


(4b) 
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and,  similarly  for  the  amplitudes 


w*  - (-)21  . (.) I+M  n 


‘M 


-M 


(5a) 


s£I3t  - (-)21  ai'> . t-)1-'1 


(5b) 
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With  these  definitions  there  holds  for  any  quantity:  A = A.  The  equations 
(3)-(5)  are  valid  even  If  the  a^1^,  a^  are  creation  and  annihilation 
operators,  i.e.,  if  the  problem  is  treated  in  the  occupation  number  formalism. 
If  they  are  simply  c-number  amplitudes,  then  there  holds 


(6a) 


which  can  be  used  together  with  eq.  (4)  and  (5)  to  define  all  amplitudes. 

r o 

Finally,  for  some  quantities,  e.g.,  =:  (-i)  Y , there  holds 


rfi]  - A[I] 


(6b) 


Such  quantities  are  called  "self  conjugate." 

With  our  conventions,  the  conjugation  (3)  becomes 

J = £(-)21 

M 

We  now  are  prepared  to  introduce  the  irreducible  density  matrices,  which 
we  shall  write  as  pjj^.  They  are  defined  in  the  reduction  of  the  direct 
product  J y into  components  of  good  angular  momentum.  We  achieve  this  by 
recoupling  (see  Fig.  1),  to  obtain 
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Figure  1 
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Here  we  have  introduced  the  irreducible  density  matrix 


by  the  definition 


^L]  • 
m L 


(8) 


Below  we  will  give  a justification  for  the  choice  of  the  coupling  order. 

Before  continuing  we  would  like  to  remark  that  particular  caution  must  be 
paid  if  the  amplitudes  a^  are  in  fact  annihilation-creation  operators.  So, 
for  example,  if  one  is  dealing  with  a many-body  system  and  if  one  writes  out 
the  operator  a^  explicitly  in  terms  of  single-particle  operators,  e.g.,  for 
a two-body  system 


4I]  - [°CjVk]j“] 


(9) 


one  has,  in  view  of 

(aB)+  = Btat 


(10) 


5 


~rn 

for  the  adjoint  operator  aj^  , i.e.,  the  creation  operator,  the  order 


=£13  _ / J+k-I,~Ck]~[j]  ,[I] 
aM  “ M Lp  a Jm 


(ID 


More  generally,  in  the  case  of  a many-body  system,  e.g.. 


A^n  - [[[W0^1^  ...]CO 


(12) 


one  has  for  the  adjoint  operator 


11  _ ^ ^ j k + £+ • • • - 1 


[...[y1-^  [^k]cf'i:']CI^][I2^]^I]  . (13) 


The  phase  factor  is  recognized  as  the  usual  phase  factor  associated  with 
the  re-ordering  of  the  coupling.  It  is  independent  of  the  details  of  the 
coupling  scheme,  i.e.,  of  the  intermediate  angular  momenta  Ilt  I2,  etc.  The 
same  phase  factor  applies  also  if  the  particles  occupy  identical  orbits,  i.e., 


rate,  usually  it  is  not  advantageous  to  employ  the  occupation  number  represen- 
tation when  one  works  with  coupled  operators  as,  for  example,  with  eq.  (12), 
but  it  is  easier  to  achieve  the  required  symmetry  of  the  wave  functions  by 
means  of  a fractional  parentage  expansion.  A detailed  discussion  of  this 
subject  is  contained  in  Ref.  5. 

We  now  can  evaluate  the  mean  value  of  an  operator  for  the  state  (4). 

Take  an  operator  of  multipolarity  K: 


Ni  - 1 

if  j = k = £,  etc.  It  then  is  (-)  J for  N particles  in  the  j shell.  At  any 


.00  oCK]  . ;,[K]o[K],[0] 


(14) 
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Then  we  have,  according  to  Fig.  2 


[a[i]  „ M a[i]] t#] 


Figure  2 


<f|n|  = l I2  K [k  K o][i  I ol  (-)1+K_I  (.)I+I‘K 
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= (-)21  [PC«JK]][°]  ^[I]|s,[K]|  JIHj  > (15) 

The  structure  of  (15)  displays  the  well-known  obvious  facts:  (i)  in  order  to 

yield  a finite  expectation  value  the  system  must  have  a density  matrix  which 
has  a non-vanishing  component  of  the  required  multipolarity;  in  other  words, 
only  one  multipolarity  of  the  density  matrix  contributes  to  the  expectation 
value,  (ii)  the  orientations  of  the  operator  (given  by  the  set  ^ and  of 
the  system  (specified  by  the  set  ) must  match;  in  other  words,  the 
expectation  value  is  proportional  to  the  quantity  , which  is  the 

generalization  to  arbitrary  K of  the  usual  scalar  product  of  vector  algebra, 
where  K = 1. 
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The  phase  (-)  in  eq.  (15)  arises  since  the  density  matrix  incorporates 
a reordering  transformation.  To  see  this  more  clearly,  consider  the  general 
case  for  pure  states,  i.e.,  the  case  of  a transition  matrix  element,  say,  for 
the  transition  Tj  + Changing  in  Fig.  2 Tj  to  y.  , I.e.,  changing  I to  L In 
the  top  two  input  lines,  we  have 


<\|il|  f,> 


= I K L 


fL  L 0“ 

"I  I O' 

K K 0 

I I 0 

Li  i o. 

.0  0 0 

[S 


scn.iK],CL]i[o]  r,.,[L]i _[k] i „ri] 


]LUJ  [fLJ|sfKJ|f 1J] 


- [iCL]jK]acn][o] 


(16) 


K+I-L  (_,  I+L-K  ^[LljtnjWjCO] 


(16a) 


The  step  from  (16)  to  (16a)  represents  the  re-ordering  shown  in  Fig.  2 
after  the  dashed  line  C. 

For  an  impure  state  one  cannot  factorize  the  amplitudes  and  one  must 

bring  the  state  amplitudes  together  in  order  to  achieve  a density  matrix.  One 

actually  has  two  choices  for  defining  the  density  matrix,  the  choice  (8), 

employed  in  Fig.  2,  and  the  opposite  order,  exchanging  a and  a.  In  that  case 

21 

the  second  crossing  in  Fig.  2 would  be  absent,  and  the  phase  (-)  of  (15) 
would  be  replaced  by  (-)  . Clearly,  the  first  choice  is  preferable: 
independently  of  the  multipolarity  of  the  operator  and  the  density  matrix,  the 
introduction  of  a density  matrix  for  a half-integer  system  yields  a minus 
sign,  i.e.,  the  value  of  the  crossing  boxes  of  Fig.  2. 

The  structure  of  eq.  (8)  suggests  the  definition  of  the  reducible  density 
matrix: 


- aCMI] 

pmm‘  * am  am' 


(17) 
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This  matrix  is  fully  contra-standard.  The  coupling  to  good  angular  momentum, 
i.e.,  the  construction  of  an  irreducible  density  matrix,  is  performed  using 
the  normal  vector  coupling  coefficients: 

P[L]=  l (Im  Im'lLM)  p , . (18) 

mm' 

A more  general  case  arises  in  the  context  of  the  evaluation  of  a 
transition  probability.  Formally  it  arises  by  writing  for  the  operator  (14) 
the  direct  product 


fl  = | T+  |f>  <f | T | . (19) 

The  operator  T can  be  a single  multipole  or,  as  in  a plane  wave,  a mixture  of 
multipoles.  This  way  (15)  then  yields  the  square  of  the  transition  matrix 
element: 


<^1  Tf  | vf>  <Yf|  T | = |Mj  2 , (20) 

which  is  directly  the  transition  probability,  of  course,  except  for  the 
kinematical  factors.  Writing 

T =:  l K [t[K^  T^  j[0:i  (21) 

K 


we  obtain  for  (21)  from  Fig.  3 
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Figure  3 


Figure  4 


10 


p ( L I 

P |HJ 
T 

■^|PJ 

[0] 


Figure  5 
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For  self-conjugate  transition  operators  one  may  omit  the  tildas  in  (22b).  In 

the  case  where  the  operator  consists  of  a single  multipole,  the  sum  over  K,K' 

collapes  into  a single  term.  Note  that  the  final  form  (22b)  still  contains  a 

recoupling  coefficient  (here,  a 6-j  type).  This  is  the  consequence  of  the 

★ 

fact  that  the  density  matrix  ties  together  the  matrix  elements  M and  M , as 

^ * 

evident  from  (8):  the  amplitude  a originates  in  M,  the  amplitude  a in  M . 

The  recoupling  coefficient  in  (22b)  unties  the  invariant  matrix  elements.  It 
thus  seems  that  for  pure  states,  the  use  of  the  density  matrix  only  introduces 
an  unnecessary  complication.  But  even  for  pure  states  the  use  of  the 
amplitude  rather  than  the  density  matrix  allows  for  no  simplification  when 
analyzing  the  final  state,  as  demonstrated  below. 

In  equation  (22)  the  information  of  the  particular  experimental  set-up, 
i.e.,  the  preparation  of  the  initial  state,  etc.  is  contained  in  the  invariant 
product  [ p^t^t^  . Of  course,  here 


[£C«  ]CLD  ..  ^L] 


(23a) 


are  the  density  matrices  characterizing  the  transition  operator. 

To  obtain  the  characteristics  of  the  complete  reaction,  e.g.,  the  angular 

correlations  in  (y,p)  or  (y»P»n)  reactions,  etc.,  one  simply  continues  the 

graph  by  inserting  an  appropriate  analysis  section.  We  illustrate  this  by  the 

case  of  the  emission  of  one  particle.  The  analysis  section  here  is  inserted 

in  the  region  shown  by  the  dashed  line  A-A  in  Fig.  3.  In  its  place  we  insert 

the  graph  shown  in  Fig.  4,  which  also  shows  the  line  A-A,  beyond  which  lies 

the  original  graph  of  Fig.  3.  The  operator  which  analyzes  the,  say,  angular 

distribution,  or  polarization,  etc.,  of  the  emitted  particle  is  [wfi];  it 

rpl 

provides  the  density  matrix  , 
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(23b) 


l JCA3jCp3 

A, A' 


We  here  have  allowed  for  the  general  case  that  the  analyzing  operator  is  given 
by  a mixture  of  multi  poles.  The  angular  momentum  P replaces  some  zeroes  in 
the  earlier  coupling  graph.  Fig.  3.  Therefore,  some  quantum  numbers  of  Fig.  3 
are  not  fixed  and  are  to  be  replaced  by  the  quantum  numbers  given  in  the  first 
parentheses;  in  general  they  must  be  summed  over  their  respective  range  of 
allowed  values.  For  example,  the  density  matrices  p and  p^  instead  of  being 
coupled  to  zero  are  coupled  to  P,  and  hence  the  density  matrix  does  not 
have  to  have  multipolarity  L;  consequently  L here  is  replaced  by  R.  In 
Fig.  4 we  also  notice  the  symbols  for  the  fractional  parentage  coefficients 
(or,  more  generally,  the  spectroscopic  factors),  which  break  the  A-body  system 
into  a one-body  plus  an  (A-l)  - body  system.  The  angular  distribution 
(and/or  polarization)  of  the  emitted  particle  (labelled  j)  is  analyzed  by  the 
operator  n,  the  orientation  of  which  is  described  by  the  density  matrix  p^. 

The  angular  distribution  etc.,  is  given  by  the  quantity  [p^  pyR^  p^]^* 
of  course,  summed  over  all  values  of  L,  R,  P,  contained  in  the  experimental 
situation.  We  shall  not  bother  to  write  out  the  algebraic  expression  of 
Fig.  4.  Of  course,  the  complete  expression  is  the  product  of  the  expressions 
of  Fig.  3 and  Fig.  4. 

We  now  continue  the  analysis  to  the  case  of  the  sequential  emission  of  a 
second  particle.  To  that  end  one  continues  the  graph  by  replacing  the  region 
B-B  of  Fig.  4 by  the  graph  of  Fig.  5.  Consequently  some  quantum  numbers  of 
Fig.  3 and  Fig.  4 must  be  generalized,  as  indicated  in  the  drawing;  in  Fig.  4 
by  the  quantum  numbers  in  the  parentheses  and  in  Fig.  3 by  the  quantum  numbers 
in  the  second  parentheses.  The  essential  point  is  that  the  wave  functions  of 
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the  ( A-l)  system  are  not  terminated  in  overlap  matrix  elements  but  further 
decomposed  into  a one-body  and  an  (A-2)-body  system.  The  one-body  part  is 
analyzed  by  the  operator  W,  and  the  (A-2)-body  system  is  overlapped  away.  The 
details  of  this  part  of  the  problem  have  been  discussed  in  Ref.  4,  and  even 
more  extensively  in  Chapter  6 of  Ref.  5.  Clearly  this  procedure  can  be 
continued  to  achieve  the  complete  description  of  reactions  resulting  in  the 
emission  of  arbitrarily  many  particles,  by  cutting  the  graph  at  C-C  and  in 
essence  repeating  the  procedure  of  Fig.  5.  In  other  words,  the  procedure  of 
going  from  n detected  reaction  products  to  n + 1 detected  reaction  products  is 
the  same  as  that  going  from  one  to  two  detected  particles,  shown  in  Fig.  5. 

The  density  matrices  here  arrive  in  the  "natural"  coupling  scheme,  i.e., 

[ [ [ [ p p-j-]  p^]  pg]  . ..J^.  This  may  not  be  the  desired  form;  for  example,  one 
may  be  interested  in  the  correlation  between  the  two  emmitted  particles.  Then 
one  would  re-couple  the  density  matrices  as  [[p  Pj]Cpa  Pg]  In  the 

extreme  case  where  one  is  interested  only  in  the  correlation,  one  would  put 
L = T = 0,  i.e.,  ignore  the  unnecessary  information.  Or,  one  wants  only  the 
total  cross  section.  Then  one  restricts  the  summation  to  the  terms  P = Q = 0; 
analogously  for  any  other  situation.  We  here  shall  not  go  into  further 
details,  as  our  aim  was  only  to  demonstrate  the  use  of  the  irreducible  density 
matrices,  and  the  ease  of  arriving  at  the  final  results  without  any  problems 
in  phasing,  in  particular  in  the  interference  terms  which  here  arise 
automatically  in  the  form  of  the  summation  over  the  angular  quantum  numbers 
allowed  by  triangularity  and  contained  in  the  experimental  setup. 
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III.  IMPURE  STATES 


If  the  system  is  not  in  a pure  state,  it  cannot  be  described  by  a wave 

function.  In  that  case,  in  essence  one  must  start  with  Fig.  2.  The  density 
riel 

matrices  p then  must  be  given  directly.  The  system  is  still  defined  by 

eq.  (7),  except  that  the  left  hand  side  of  (7)  has  only  symbolic  meaning. 

This  way  we  arrive  at  the  following  rule  of  working  with  density  matrices 

for  impure  systems.  One  may  begin  with  the  output  side  of  Fig.  2,  i.e.,  with 

the  right-hand  side  of  eq.  (16).  However,  in  case  of  doubt,  make  a graph  for 

the  case  of  pure  states.  Then  introduce  the  permutation  in  the  invariant 

triple  product  of  the  amplitudes  as  in  Fig.  2,  and  replace  the  "pure"  density 

matrix  [a^a^ ^ by  the  density  matrix  describing  the  actual  system. 

FkI 

That  means,  that  from  now  on  for  impure  systems  pL  J must  be  treated  as  a 

single  tensor.  Otherwise  all  usual  rules  of  writing  down  the  algebraic 

expressions  for  the  recoupling  diagram  remain  unchanged.  In  particular,  the 

central  recoupling  graph  rule  holds:  write  down  the  product  of  all  symbols 

appearing  in  the  graph,  including  all  crossing  phases.  Precisely  for  this 
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reason  it  is  inadvisable  to  incorporate  the  phase  (-)  of  eq.  (15)  into  the 
definition  of  the  density  matrix,  even  though  the  eq.  (15)  then  would  look 
simpler  in  that  it  would  have  no  explicit  phase.  However,  then  the  central 
rule  of  using  recoupling  graphs  would  be  broken. 

As  an  example  for  an  impure  system  consider  a cryogenic  alignment 
experiment.  A nucleus  interacts  at  temperature  T with  the  crystalline  field. 
Let  us  assume  that  this  field  is  given  by  a mixture  of  multi  poles.  The 
Hamiltonian  then  can  be  written 


We  first  determine  the  eigenstates  of  the  Hamiltonian,  which  of  necessity  are 
pure  states.  Thus,  writing 


H*n 


E ¥ 
n n 


(25) 


and,  if  I is  the  spin  of  the  nucleus. 


Yn  = i [s^  J-1-^  j'-0-' 


(26) 


Fi gure.6 


we  have  (see  Fig.  6) 


l K I 

MK 


K K O' 
I I 0 
I I 0 


"I  I 0" 
0 I I 
.1  0 I. 


X (IMI  -Mj  00)  [hCKlSCnI]]CMn 


■ I En<IMI  Sn” 


(27) 
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Equating  the  coefficients  of  and  evaluating  the  recoupling  coefficients 
we  obtain 


l [/I]|H[I°|^I]]  (KkIM'|IM)  h[K]S^?  = E S™ 
KkM'  * 


(28) 


which  is  the  eigenvalue  equation  for  y . Using 


cCl]  _ , ,I-M  S[I]* 
^nM  ‘ ^n-M 


one  can  construct  the  density  matrix  for  the  pure  states 


anMM' 


cCl]  c{I] 

^nM  anM'  * 


(29) 


The  density  matrix  for  the  target  nucleus  at  equilibrium  with  the  target 
crystal  at  temperature  T = 1/3  is  given  by 

- 3En  _ 

^MM'  = 2-1  ^ e anMM' 
n 

-SEn 

z = l e . (30) 

n 

One  sees  that,  indeed,  p^,  cannot  be  factorized.  It  can  only  be  written  as  a 
sum  of  factorizable  matrices. 

The  irreducible  density  matrices  are  again  constructed  by  eq.  (18).  They 
are  handled  exactly  as  the  irreducible  density  matrices  for  pure  states. 
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Of  course,  one  encounters  cases  where  the  system  is  not  in  a state  of 
definite  angular  momentum.  The  most  prominent  of  such  cases  is  a scattering 
state  where  the  incoming  particle  is  a plane  wave.  No  difficulties  arise  in 
this  case,  as  has  been  already  discussed  above.  The  independent  summation 
over  the  different  multipoles  ( K , K' , A,  A',  etc.)  yields  directly  the  inter- 
ferences between  the  different  multi  poles;  there  is  no  difference  between  the 
cases  of  the  pure  and  the  impure  states,  when  the  expressions  are  written  in 
terms  of  the  density  matrices. 
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IV.  RELATION  TO  THE  USUAL  DEFINITIONS 


The  density  matrix  usually  is  defined  by 


. ■ a'1*  a<!>*  . (-)21  a'1’  & 


mm  m m 


m nr 


(31) 


and  the  state  is  given  by 


p , /I]  w 

Kmm ' rm  ym 


(32) 


The  mean  value  of  a component  of  an  operator  then  is 


<!5kK]> 


= Tr  p n 


= l p , 

L - Mmnr  V 1 k 1 ym 


mm 


(33) 


Here  both  the  matrix  elements  and  the  density  matrix  are  mixed  quantities  in 
both  gender  and  angular  momentum.  We  have. 


/ \ I-m  - 

pmm'  " p-mrn' 


(34) 


and  thus  the  reduction  in  angular  momentum  is 


P„L]  = l (I-m  Im'|LM)  (-l1""  p^, 


(35) 
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which  is  not  the  familiar  form  for  coupling  angular  momenta.  When  using 
density  matrices,  it  is  therefore  advisable  to  work  exclusively  with  p.  As 
demonstrated  above  with  the  example  of  the  polarized  target  it  is  always 
possible  to  obtain  directly  the  fully  contrastandard  p.  There  thus  is  no  need 
to  introduce  at  all  the  mixed  gender  p.  At  any  rate,  if  one  has  p available 
one  can  construct  p by  means  of  (34). 
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APPENDIX 


For  the  convenience  of  the  reader  we  list  here  some  symbols.  More 
details  are  given  in  Chapter  4 of  Ref.  5. 


1.  Definition  of  invariant  matrix  element 


r J!j  ^]+  y[K]  Ij  2]  \ J 2 f ^1  ^ J 2 ) 

J mL  k m2  ['-ml  k 


| y[K]  j ^J]j  21 


J.  (Al) 


2.  Relation  between  the  invariant  matrix  element  and  the  reduced  matrix 
el ement 


< i]  j |j[K]  1 1 ^j2]>  = (_)J’i+K-j2  [/J  i.]  |t^K]  | ^J2]]  . ( A2) 


3. 


Special  case 


= l: 


overl  ap 


[/I]|l|  Z1-*]  =:  [Z1-*!  = I 


(A3) 


4.  Amplitude  overlap 


[anvi]j  - 


1/1 


( A4) 


5.  9j  - type  recoupling  coefficient 


'a  b c' 
d e f 
_g  h i. 


( A5) 
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